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Abstract 



Given a set of points that sample a shape, the Rips complex of the data points is often used in 

machine-learning to provide an approximation of the shape easily- computed. It has been proved 

^N recently that the Rips complex captures the homotopy type of the shape assuming the vertices of 

the complex meet some mild sampling conditions. Unfortunately, the Rips complex is generally 

I— I high-dimensional. To remedy this problem, it is tempting to simplify it through a sequence of 

rn collapses. Ideally, we would like to end up v^rith a triangulation of the shape. Experiments suggest 

f) that, as we simplify the complex by iteratively collapsing faces, it should indeed be possible to 

• avoid entering a dead end such as the famous Bing's house with two rooms. This paper provides 

Q a theoretical justification for this empirical observation. 

'— ' We demonstrate that the Rips complex of a point-cloud (for a well-chosen scale parameter) 

I can always be turned into a simplicial complex homeomorphic to the shape by a sequence of 

^ collapses, assuming the shape is nicely triangulable and well-sampled (two concepts we will ex- 

f^ plain in the paper) . To establish our result, we rely on a recent work which gives conditions under 

OO which the Rips complex can be converted into a Cech complex by a sequence of collapses. We 

proceed in two phases. Starting from the Cech complex, we first produce a sequence of collapses 

that arrives to the Cech complex, restricted by the shape. We then apply a sequence of collapses 

'^ that transforms the result into the nerve of some "robust" covering of the shape. 
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1 Introduction 

This paper studies tlie problem of converting a Cecil complex whose vertices sample a shape into a 
triangulation of that shape using collapses. 

Imagine we are given a set of points that sample a shape and we want to build an approximation 
of the shape from the sample points. A practical approach, often used in machine learning, consists 
in outputtingthe Rips complex of the points [8, 13, 12]. Formally, the i?/p5compZex of a set of points P 
at scale a is the abstract simplicial complex whose simplices are subsets of points in P with diameter 
at most 2a. The Rips complex is an example of a flag complex — the maximal simplicial complex 
with a given 1 -skeleton. As such, it enjoys the property to be completely determined by its 1 -skeleton 
which therefore offers a compact form of storage easy to compute. Recently, precise sampling con- 
ditions have been formulated under which the Rips complex reproduces the homotopy type of the 
shape [9, 3,5]. Unfortunately, the dimension of the Rips complex can be very large. This suggests a 
two-phase algorithm for shape reconstruction. The first phase builds the Rips complex of the data 
points, thus producing an object with the right homotopy type. The second phase simplifies the Rips 
complex through a sequence of collapses. Ideally, we would like to end up with a simplicial complex 
homeomorphic to the underlying shape. 

Yet it is not at all obvious that the Rips complex whose vertices sample a shape contains a subcom- 
plex homeomorphic to that shape. Even if such a subcomplex exists, is there a sequence of collapses 
that leads to it? Certainly if we want to say anything at all, the geometry of the complex will have to 
play a key role. As evidence for this, consider a simplicial complex whose vertex set is a noisy point- 
cloud that samples a 0-dimensional manifold and suppose the complex is composed of a union of 
Bing's houses with two rooms, one for each connected component in the manifold. Each Bing's house 
is a 2-dimensional simplicial complex which is contractible but not collapsible. Thus, the complex 
carries the homotopy type of the 0-dimensional manifold but is not collapsible. Fortunately, it seems 
that such bad things do not happen in practice, when we start with the Rips complex of a set of points 
that samples "sufficiently well" a "nice enough" space in R'^ . The primary aim of the present work is 
to understand why. For this, we will focus on the Cech complex, a closely related construction. For- 
mally, the Cech complex of a point set P at scale a consists of all simplices spanned by points in P 
that fit in a ball of radius a. In [5] , it was proved that the Rips complex can be reduced to the Cech 
complex by a sequence of collapses, assuming some sampling conditions are met. 

In this paper, we give some mild conditions under which there is a sequence of collapses that 
converts the Cech complex (and therefore also the Rips complex) into a simplicial complex home- 
omorphic to the shape. Our result assumes the shape to be nicely triangulable, a concept we will 
explain later in the paper. Perhaps unfortunately, the proof of existence is not very constructive: it 
starts by sweeping space with offsets of the shape — which are unknown — and builds a sequence of 
complexes which have no reason to remain close to flag complexes and therefore cannot benefit from 
the data structure developed in [4] . Nonetheless, even if results presented here do not give yet any 
practical algorithm, we believe that they provide a better understanding as to why the Cech complex 
(and therefore the Rips complex) can be simplified by collapses and how this ability is connected to 
the underlying metric structure of the space. 

We now list the principal results of the paper, materialized as brown arrows in Figure 1. In Sec- 
tion 3, we study how the reach of a shape is modified when intersected with a (possibly infinite) 
collection of small balls. Based on this knowledge, we introduce the Cech complex restricted by the 
5/jflpe A in Section 4 and give conditions under which it is homotopy equivalent to A (Theorem 1). In 
Section 5, we find conditions under which there is a sequence of collapses that goes from the Cech 



complex to the restricted Cech complex (Theorem 2). Combined with Theorem 1, this gives an alter- 
native proof to a result [15] recalled in Section 2 (Lemma 2). In Section 6, we define a-robust coverings 
and give conditions under which the restricted Cech complex can be transformed into the nerve of 
an a-robust covering (Theorem 3). In Section 7, we define and study nicely triangulable spaces. Such 
spaces enjoy the property of having triangulations that can be expressed as the nerve of a-robust 
coverings for a large range of a. Finally, we exhibit examples of such spaces. Section 8 concludes the 
paper. 

Rips(i^a) ^^^-^w-^ws^ Cech(j;a) ---^^•~^~~'^~-> Cech^(j;a) -^-^•^•^•■-^'-^^^ Nerve of an a-robust covering of A 

Lemma 1 
(Nerve Lemma) 

Lemma 2 [15] 

"^ — >A 




Figure 1: Logical structure of our results. Brown arrows represent new results. The arrow ^^ stands for "deformation 
retracts to". The arrow •"» stands for "can be transformed by a sequence of collapses into". The symbol "^" means 
"homotopy equivalent to" and "i^" means "homeomorphic to". 



2 Background 

First let us explain some of our terms and introduce the necessary background. M^ denotes the d- 
dimensional Euclidean space. ||x — y|| is the Euclidean distance between two points x and y of M^. 
The closed ball with center x and radius r is denoted by B[x, r) and its interior by B°[x, r). Given a 
subset XczR'', the a-offset of X is X®« = [j^^^ B{x, a). The Hausdorff distance dnlX, Y) between the 
two compact sets X and Y of W^ is the smallest real number e >0 such that X c Y®'^ and Y c X®'^. 
We write d[x, Y) = infygy ||x — y\\ for the distance between point x e M^ and the set Y czM.'^ and 
d{X, F) = inf;t:6xinfyey||x-y|| for the distance between the two sets X c R^ and FcR'^. 

A convenient way to build a simplicial complex is to consider the nerve of a collection of sets. 
Specifically, let 'i^ = {Cp | p e P} be a 1 -parameter family of sets indexed by p e P. The nerve of 
the family is the abstract simplicial complex that consists of all non-empty subcoUections whose sets 
have a non-empty common intersection, Nrv*)^ = {cr c P | cr 7^ and npecr ^p 7^ ^^ ^^ ^^^^ paper, we 
shall consider nerves of coverings of a shape A. We recall that a covering of A is a collection "ia ={Cp\ 
p e P} of subsets of A so that A = (JpeP ^p- ^^ ^^ ^ closed {resp. compact] covering if each set in "iS is 
closed {resp. compact) . It is a finite covering if "^ is finite. The Nerve Lemma gives a condition under 
which the nerve of a covering of a shape shares the topology of the shape. It has several versions [6] 
and we shall use the following form: 

Lemma 1 (Nerve Lemma). Consider a compact set A c M^. Let "iS" = {Cp \ p ^ P] be a finite closed 
covering of A. If for every 7^ cr c P, the intersection Hzeo- ^^ is either empty or contractible, then the 
underlying space o/NrvS^ is homotopy equivalent to A. 

Hereafter, we shall omit the phrase "the underlying space of" and write X^^Yto say that X is ho- 
motopy equivalent to Y. Given a finite set of points P e R^ and a real number a > 0, the Cech complex 
of P at scale a can be defined as Cech(P a) = Nrv{B{p, a) | p e P}. With this definition and the Nerve 
Lemma, it is clear that Cech(P a) ~ P®" ; see the black vertical arrow in Figure 1 . Several recent results 



have expressed conditions underwhich P®" recovers the homotopy type of the shaped [15, 11, 10, 5]. 
Intuitively, the data points P must sample the shape A sufficiently densely and accurately. One of the 
simplest ways to measure the quality of the sample is to use the reach of the shape. Given a compact 
subset AoiR'^, recall that the medial axis Ma of A is the set of points in M^ which have at least two 
closest points in A. The reach oiA is the infimum of distances between points in A and points in Ma, 
Reach(A) = infaeA.me^^^ ll*^ ^ '^ll- We have (horizontal black arrow in Figure 1): 

Lemma 2 ([15]). Let A andP be two compact subsets ofR'^ . Suppose there exists a real number e such 
thatdniA, P) < s < (3 - ^8) Reach (A). Then, P®" deformation retracts to A for a = [2 + V2)e. 

Combining Lemma 1 and Lemma 2, we thus get conditions under which Cech(P a) :^ A. The next 
two sections will provide an alternative proof of this result along the way. 

Before starting the paper, let us recall what a collapse is. Consider a simplicial complex K and 
CTmin £ K. Let A be the collection of simplices in K having amin as a face. Provided that there is a 
unique inclusion-maximal element CTmax 7^ CTmin in A, it is well-known that K deformation retracts to 
^\ A and the operation that removes A is called a collapse. 

3 Reach of spaces restricted by small balls 

In this section, we consider a subset A czW^ such that Reach [A] > and prove that as we intersect 
A with balls of radius a < Reach (A) the reach of the intersection can only get bigger; see Figure 2. 
More precisely, let A be a compact subset of M^ and a c R^. Write dS[a,a) = Hzeo- ^i^'^^ ^^^ the 
common intersection of balls with radius a centered at a and assume that A n ^(cr, a) 7^ 0. In this 
section, we establish that Reach (A) < Reach (A n ^(cr, a)) first when a is reduced to a single point z 
(Lemma 3), then when a is finite (Corollary 1) and finally when cr is a compact subset of M^ (Lemma 
5). Although the first generalization is all we need for establishing Theorem 1 in Section 4, the second 
generalization will turn out to be useful later on in the paper. 
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Reach (A nC) 



Figure 2: Left: Medial axis ^^ of a shape A. If we intersect A with a ball B whose radius is smaller than the reach 
of A, the reach of the intersection AUB can only get bigger. Right. This property does not hold if we replace the 
ball B by another set, even with infinite reach such as the cube C. 



Lemma 3. Let A c M^ be a compact set and B{z, a) a closed ball with center z and radius a. IfO <a< 
Reach (A) and An B{z, a) 7^ then Reach (A) < Reach (A n B{z,a)). 



Proof. See Figure 3 on the left. Set a^ = a and let A' = An B[z,az). Assume, by contradiction, 
that Reach(A') < Reach(A) and consider a point y in the medial axis of A' such that d[y,A') = Uy < 
Reach [A). Introduce A" = A'r\ B{y, ay ) and denote by c and p the center and the radius of the small- 
est ball enclosing A". Because A" is contained in both B[z,az) and B{y, ay), the radius of the smallest 
ball enclosing A" satisfies p < min{az,ay}. Because A" c A' c A we get d{c,A) < d[c,A') < d[c,A") < 
p <ay < Reach (A) and therefore c has a unique closest point tzaIc) in A. Take r to be any real number 
such that max{az,ay} < r < Reach(A). We claim that r — y r^ — p^ < d{c,A); see Figure 3 for a geo- 
metric interpretation of the quantity r - \l r'^- p"^. Indeed, for every m e {z,y\, since A" c B[rn,am) 
and am > p. Lemma 12 (i) implies that the following inclusion holds: 



B(c,am--^a^-p2) c B(m,am). 
Since the map r — > r — y/r^ — p2 is strictly decreasing in [p, +oo), we get that 

B(c,r-Vr2-p2) c B°(z,a^)nB°(y,a3,). 

By definition, points of A in B°(z,az) belong to A' and B°(y,ay) contains no points of A'. It follows 
that B(c, r — \J r'^ — p^) contains no points of A and r — \j r"^ — p^ < d{c,A) = \\c — 7r^(c)|| as claimed. 
Let us consider the point x = 7ryi(c) + ,,'^ J c — 7ryi(c)); see Figure 3 on the right. By construction 
||x - 7ryi(c)|| = r < Reach(A) and therefore x has a unique closest point nA[x) = UaIc) in A. Since 
B(c, £i(c. A)) c B(x, ti(x, A)), we deduce that 

B(c,r--^r2-p2) c B°(x,r). 

Applying Lemma 12 (ii) we get that A" n B°(x, r) 7^ and therefore B(x, r) contains points of A in its 
interior. But this contradicts d[x,A) = r. O 





Figure 3: Notation for the proof of Lemma 3. The quantity r- \/r^- p^ represents the height of a spherical cap 
whose base has radius p and which lies on a sphere with radius r. 



Corollary 1. Consider a compact set A c R^ and a finite set a c E^. IfO <a< Reach(A) andAn 
dB{a, a) 7^ then Reach(A) < Reach(A n dB{a, a)). 

Proof. By induction over the size of cr. D 

The following lemma is a milestone for the proof of Lemma 5. 

Lemma 4. Let {An)neN be a sequence of non-empty compact subsets ofR^ decreasing with respect of 
the inclusion order If there exists a real number r such that < r < Reach(A„) for all n e N, then 
r< Reach (n„gNA„). 

Proof. Letting A = H^sn^"' *^ ^^^^ show that the Hausdorff distance dniAniA) tends to as n ^ 
+00. For f > 0, introduce the set If = {x e M^ | d[x,A)>s} and notice that P„gj^(LfnA„) = I(;n^ = 0. 
Since the sequence (L^ nA„)„gN is decreasing and consists of compact sets, the only possibility is that 
LgDAi = for some i e N. Equivalently, dHiAi,A) < e which proves the convergence of A„ to A under 
Hausdorff distance. 

We now claim the following. For any positive real numbers 5 and r' in the open interval (0, r) and 
any point z' whose distance to A is r', we can find a point c such that 

AnB{z',r')(zB[c,5). 

The claim implies that A n B(z', r') is a singleton and this proves the lemma. We prove the claim by 
providing a construction for point c; see Figure 4, left. For this, take a real number e such that: 



0<e<[r-r') 



f \ 

1 

- -1 



(1) 



l-(f) ] 



Since dniAntA) — > as n — > +cxd, there is some i e N such that dniAiyA) < e. By hypothesis, r < 
Reach(A;). Since A c Ai, one has d[z',Ai) < d{z',A) = r' < r < Reach(A/) and therefore z' has a 
unique closest point at in At. Let us define the point z by 

Since r < Reach(A/), the open ball B°(z, r) is disjoint from A, and since A c At we obtain that A n 
B°(z, r) = from which we deduce 

Ar\B[z',r')<zB[z',r')\B''[z,r). 

Let c be the centre of the (li — 2)-sphere S which is the intersection of the two [d — l)-spheres d B[z', r') 
and d B[z, r). Observe that c lies on line zz'. Take a point x e S and define 9 = Zczx. It follows from 
the above inclusion that A n B[z', r') c B(c, r sin 0). Hence, to prove the claim it suffices to establish 
that r sin 6 <5. Define y as the orthogonal projection of z' on line xz. We have 

cos9 = -— -. 

\\z'-z\\ 

Observe that \\y — z\\ = \\x — z\\ — ||x — y|| > ||x — z\\ — \\x — z'\\ = r — r'. Since dniAyAt) < s, one has 
d{z,A) < d{z,Ai) + e = r + e and since d{z',A) = r' one has ||z' - z|| < r -\-e — r' . This gives using 
Equation (1) 

cos6)> '""'" , >t/l-f-l (2) 

r-\-e — r y \r ) 

from which we deduce immediately that r sin 6 <5 as required. D 





Figure 4: Notation for the proofs of Lemma 4 (left) and Lemma 5 (right). 

Lemma 5. Let A and a ^% be two compact sets ofW^. IfO < a < Reach (A) and An d^[cr,a)^% then 
Reach(A) < Reach(A n m[cT, a)). 

Proof. Consider a sequence {z/l/g^ which is dense in a. Corollary 1 implies that for all n > 0, we 
have Reach f An P"^jB(z/,a)j > a. Applying Lemma 4 we get that ReachfAnf^jgj^B(z/,a)) > a. 
We claim that AnP|j.gji^B(z,-,G) = Anf^^g^B(z,a). One direction is trivial. If a e Anf^^g^B(z,a), 
then a e Ar\[~\-^^V>[zi,a). For the other direction, take a e Anf^j.gj^B(Z(,a) and let us prove that 
Vz e (T, \\z — a\\<a. Assume, by contradiction, that for some z^cr, one has \\z — fl|| > a. Then, there 
is i eN such that ||z— z/|| < ||z — a|| — a, yielding||z/ — fl|| > ||z — a|| — ||z — z;|| > a, which is impossible. 
We have just shown that a e A n nzeo- B(z, a). D 

4 The restricted Cech complex 

Given a subset A c M"^, a finite point set P and a real number a > 0, let us define the Cech complex 
of P at scale a, Cechyi( JJ a), restricted by A as the set of simplices spanned by points in P that fit in a 
ball of radius a whose center belongs to A. Equivalently, Cech^(i^ a) = Nrv{A n B(p, a) | p e P}. In this 
section, we give conditions under which A and Cech^( J! a) are homotopy equivalent. Precisely: 

Theorem 1. Let A c M'' be a compact set, P c M^ a finite point set and a a real number such that 
< a < Reach (A) and A c P®". Then, Cec\\A{P, «) and A have the same homotopy type. 

Proof. Since A c P®", clearly A = [Jp^piA n B[p, a)). By Lemma 7 below, for all 7^ cr c P, the inter- 
section P|^g^(An B(z, a)) is either empty or contractible. We conclude by applying the Nerve Lemma 
to the collection {A n B{p, a)\peP}. D 

Recall that d3{a,a) = O^^^j B[z,a) is the common intersection of balls with radius a centered at 
a. A key argument in the proof of Theorem 1 is the fact that An ^(cr, a) is either empty or contractible, 
whenever < a < Reach (A). To establish this fact, we will need the following lemma. Let Radius (X) 
designate the radius of the smallest ball enclosing the compact set X. 
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Lemma 6. IfX c M^ is a non-empty compact set with Radius [X) < Reach [X), then X is contractible. 

Proof. We recall that for every point m such that d[m,X) < Reach (X) there exists a unique point 
of X closest to m, which we denote by Jixim). Furthermore, we know from [14, page 435] that for 
< r < Reach(X) the projection map tix onto X is I ^^l^^tL-r )-Lipschitz for points at distance less 
than r fromX. Denote by c the center of X; see Figure 4, right. If x e X and t e [0, 1], one has 

di{l-t)x + tc,X) < Wil - t]x + t c - x\\ < \\c-x\\ < Radius(X) < Reach(X). 

Thus, the map H : [0, 1] x X ^ X defined by H[t,x) = 7ix{{l — t)x + tc) is Lipschitz and defines a 
deformation retraction of X onto {7Zx{c)}. D 

We deduce immediately the following lemma. Besides being useful for proving Theorem 1, it will 
turn out to be a key tool in Section 6. 

Lemma?. Let A be a compact set ofR'^ and a a real number suchO < a <Reach{A). For all non-empty 
compact subsets a czR^, the intersection An 93[a,a) is either empty or contractible. 

Proof. Suppose A n d3[cT, a) 7^ 0. By Lemma 5, 

Radius(An^(o-,a)) < a < Reach(A) < Reach(An^(cr,a)). 
By Lemma 6, An ^{cr, a) is contractible. D 

This lemma can be seen as a variant of Lemma 7 in [1], Proposition 12 in [7] and the local reach 
lemma in [2] which all say that HA is a A:-manifold that intersects a ball B with radius a < Reach(A), 
then A n 5 is a topological fc-ball. 

5 Restricting the Cech complex by collapses 

In this section, we state a condition under which there exists a sequence of collapses that transforms 
Cech(i^ a) into its restricted version Cech^( JJ a). First we establish several facts about a collection of 
balls whose centers are close to the shape A and which have a non-empty common intersection that 
does not intersect A. These facts are formalized in the following lemma. As before, we let dS[cT,a) 
denote the common intersection of balls with radius a centered at a and by convention, we set 
d[A,@) = +00. Hence, when we write that d{A, 9S[a,a)) = t for some r e M, this implies implicitly 
that ^(cr, a) 7^0. 

Lemma 8. Let A c M'' be a compact set, a czW^ a finite set and a>0 such that d[A, dS[a,a))= t for 
some f e M. IfO <t< Reach [A) — a, we have the following properties (see Figure 5, left): 

• Thereexistsa uniquepointx e ^{cr,a) whose distance to A is t ; 

• Thesetao = {p ecr | xe d B[p,a)} is non-empty ; 

• d[A,dS{ao,a))=t. 




Clpeao^M^ 



/ P' 




HAz) 



Figure 5: Notation for the proofs of Lemma 8 and Theorem 2. Black dots belong to a and the ball B[z,a) is 
instrumental in proving Lemma 8. We show that ,SS(cro,a)c B[z,a) (on the left) using the fact that npeo- ^x(p)c 
Hji(z) (on the right), where H^[m) designates the half-space containing m whose boundary passes through x and is 
orthogonal to mx. 



Proof. Since A and 0§{a,a) are both compact sets, there exists at least one pair of points {a,x) e 
A X 0S[a, a) such that ||fl — x|| = f ; see Figure 5. By definition, x belongs to B[p, a) for all p E cr. Since 
X lies on the boundary of dS{(j, a), it lies on the boundary of B(p, a) for at least one p € cr, showing 
that ctq 7^ 0. Since Reach(^®'^) = Reach(A) — t>a, the point z=x + a ,,^~" has a unique closest point 
X in^®*^ and the distance d[A,B[z,a))= t isrealizedby the pair of points (a, x). To prove that x is the 
unique point of S^{(j, a) whose distance to A is t, it suffices to show that S^{a, a) c B(z, a). Actually, 
we will prove a stronger result, namely that ^(ctq, a) c B(z, a) which will also imply the third item of 
the lemma, that is, d{A, d3[(Jo, a))=t. 

Let us associate to every point m e M^ the half-space Hx{m) = {u e M"^, (m —x,u} > 0} and 
observe that for all m e ctq U [z] we have \\m — x\\ = a; see Figure 5, right. It is easy to check that 
since x realizes the smallest distance between points of 9B[a, a) and a, the following inclusion must 
hold Clpeao^^^P"^ '- ^x[z). This implies that for all y e Opeao^^^P"^' ^^^ distance between y and 
the boundary of Hxiz) is always larger than or equal to the distance between y and the boundary of 
Hxip) for some p^cjq. Formally, this means that for all u e M"^ and all 5 > 0, the following implication 
holds: 

-x,u)>5 =^ {z — x,u}>5 



2 



min {p 

pecro 

in the above implication and noting that for all m e aoU{z}, the following inequality 



Plugging 5: 

2{m —x,u)> ||u|p can be rewritten as \\[m —x) — u\\^< a^, we get that for all u 

implication holds: 



maxllp 

pecTo 



-X||2 <g2 



-xl|2 <a^. 



Equivalently, dS[ao, a) c B(z, a), as required. 



the following 



n 



Our second theorem (horizontal brown arrow in Figure 1) can be seen as a combinatorial version 
of Lemma 2. 
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Theorem 2. Lets > 0, a > 0, and r > 0. Consider a compact set A c M"^ li'/t/j Reach (^) > r. LetP <z W^ 
be a finite set such that dniA, P)<e. There exists a sequence of collapses fi-om Cech( J! a) to CechAiP, a) 
whenever e, a and r satisfy the following two conditions: 

(i) V2a <r-e; 

(ii) r — \l{r — sY — a^ <a — e. 

In particular, for e <{?> — ^/Q)r and a = (2 + ^/2)s, conditions (i) and (ii) are fulfilled. 

Proof. Letting ;S = r — ^{r — eY — a^, we observe that condition (i) is equivalent to ^ < r— a assuming 
a > and condition (ii) is equivalent to P < a — e. For t >0,we define the simplicial complex Kt = 
Nrv{A®f n B[p, a)\ p ^P]. Notice that Kq = Cech^(i^ a) and K+oo = Cech(i^ a). We prove the theorem 
in two stages: 

(a) First, we prove that Kt does not change as t decreases continuously from +00 to ;8. In other 
words, Kt = Cech(j; t) for all t > p. Note that this is equivalent to proving that for all non-empty 
subsets cr c P and all t > yS, 

f] B(p,a)7^0 ^^A®'nf] B(/9,a)7^0. 

pea pea 

One direction is trivial: if a point belongs to the intersection on the right, then it belongs to the in- 
tersection on the left. If the intersection on the left is non-empty, then it contains the center z of the 
smallest ball enclosing a and Radius (cr) <a< V2a <r — s. Let us apply Lemma 14 in [3] which says 
that under these conditions, the convex hull of a and therefore z cannot be too far away from A, pre- 
cisely that z e Conv(cr) c A®K Hence, z also belongs to the intersection on the right which therefore 
is non-empty. 

(b) Second, we prove that as t decreases continuously from p to 0, the only changes that may occur 
in Kt are collapses. Specifically, let Af be the set of simplices that disappear at time t, that is 

Af = {o- c P I d{A, dSia, a)) = t}. 

Suppose At 7^ and let us prove that the deletion of simplices Af from ^f is a collapse for all t e (0, p] . 
Generically, we may assume that the set of simplices Af has a unique inclusion-minimal element 
Cmin- We will explain in the appendix how to get rid of this genericity assumption. Since < t < 
P < r — a < Reach(A) — a, Lemma 8 implies that there exists a unique point x E S^{(Jmm, «) whose 
distance to A is t; see Figure 5, left. It is easy to see that Af has a unique inclusion-maximal element 
Cmax = {p e P I X e B(p, a)]. Hence, Af consists of all cofaces of CTmin and these cofaces are faces of 
o'max- To prove that removing Af from Kt is a collapse, it suffices to establish that cr^i^ 7^ CTmax- By 
Lemma 8, we know that aQ = {p^ cr min I x e 5 B[p, a)} is non-empty and belongs to Af . By the choice 
of CTmin as the minimal element of Af , we have cr min c cro and therefore x lies on the boundary of 
B{p, a) for all p e cr min- Since d{x,A) = t < p < Reach(A), there exists a unique point a e A such that 
||a — x|| = f. Because dniAyP) < e, we know that there exists a point (7 e P such that Hq' — a|| < f. Since 
||^-x|| < ||^-fl|| + ||a-x|| < e + f < e + ^ < o;, we get that x lies in the interior of 5(^, a). Therefore, 

q belongs to CTmax but not to CTmin- Hence, cr min 7^ cr max- □ 
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6 Collapsing the restricted Cech complex 

In this section, we find conditions under which there is a sequence of collapses transforming the 
restricted Cech complex Cech^(i^ a) into the nerve of an a-robust covering of ^. We start by defining 
a-robust coverings and state our result in Section 6.1. Our proof technique consists in introducing a 
2-parameter family of compact sets & = {Dp[t) \[p,t)^Px [0, 1] } and monitoring its evolution as the 
parameter t increases continuously from to 1. In Section 6.2, we give some general conditions on 
® that guarantee that the simplicial complex K[t) = Nrv{Dp(f ) | p e P} only undergoes collapses as t 
increases from to 1. We believe that these conditions are sufficiently general to be applied to other 
situations and therefore are interesting in their own right. Armed with this tool, we establish our third 
result in Section 6.3, that is, we find a family of compact sets ® which enjoys the properties required 
in Section 6.2 and such that K{0) = Cech^(i^a) and K{1) is isomorphic to the nerve of an a-robust 
covering of A. 

6. 1 Towards the nerve of a-robust coverings 

To state and prove our third theorem, we need some definitions. Given a subset X c M"^ , we call 
the intersection of all balls of radius a containing X the a-hull ofX and denote it by Hulla(X). By 
construction, Hulla(X) is convex and Hull+oo(X) is the convex hull of X. Setting ClencherSa(X) = {z e 
E^ I X c 5(z, a)}, we have 

Hull„(X)= f] Biz, a). 

zeClencherSa(X) 

Notice that Clenchers„(X) is also convex. Indeed, if two balls B{zi,a) and B{z2,a) contain X, then 
any ball B(AiZi + A2Z2, a) with Ai + A2 = 1, Ai > and A2 > also contains X. Furthermore, if X is 
compact, so is ClencherSa(X). 

Definition 1 (a-robust coverings). A covering"^ = {Cp \ v ^V] of A is a-robust if(l) each set in S^ can 
be enclosed in an open ball with radius a; (2)NTv'iS = Nrv{AnHull«(C[,) \ v ^V}. 

Of course, one may wonder if a-robust coverings of a shape A often arise in practice. Section 7 
will address this issue. For now we focus on establishing properties of a-robust coverings. 

Lemma 9. If (a is a finite compact a-robust covering of A andO <a< Reach[A), thenNrv'iS ~A. 

Proof We apply the Nerve Lemma to the coUecfion [A n Hulla(Cy) \ v eV}. Clearly, A = {J^f^y^A n 
Hulla(Cy)). By Lemma 7, for all 7^ cr c 1^, the intersection A n f^j^g^Hulla(Ci,) is either empty or 
contractible. D 

Combining the above lemma and Theorem 1 we thus get that NrvS^ ~ Cech^(J)a) for all finite 
compact a-robust coverings 'i = {d, \ v ^ V} oi A with < a < Reach (A). Next theorem strengthens 
this result and states mild conditions on P and V under which there exists a sequence of collapses 
transforming Cechyi(J! a) into a simplicial complex isomorphic to Nrv*!^. 

Tlieorem 3. Let A be a compact set ofW^ and a a real number such that < a < Reach(A). Let 
^ = \Cy \ V ^ V} be a compact a-robust covering of A. Let P be a finite point set and suppose there 
exists an injective map f -.V ^> P such that Cy c B°(/(i'), a) for all v ^V. Then, there exists a sequence 
of collapses from Cech^(i^ a) to fiNrv'^) = [fia) \ a e NrvS^}. 

The proof is given in Section 6.3. 
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6.2 Evolving families of compact sets 

In this section, we present a tool that will be useful in Section 6.3 for establishing Theorem 3. Consider 
a covering of a topological space A and suppose this covering evolves over time. We state conditions 
under which the evolution of the nerve of this covering only undergoes collapses. Conditions are 
formulated in a very general setting. We only require A to be a compact topological space: no need 
for A to possess an Euclidean structure nor a reach. We do not even need to endow A with a metric 
structure. Before stating our lemma, let us introduce one additional piece of notation. Given a finite 
set a and a map 5 : cr — > [0, 1), we write s' >~ s to designate a map s' -.a ^ [0, 1] such that s'[p) > s[p) 
for all p e a. We will say that the map 5 is constant if s[p) = s[q) for all [p, q)^a^. 

Lemma 10. Let A be a compact topological space and P a finite set. Consider a 2-parameter family 
of compact subsets of A, & = {Dp{t) \{p,t)^ P x [0, 1] }, and suppose the following four properties are 
satisfied: 

(a) ForallO< t <t' <l andallp^P, we have Dp{t') (Z Dp{t)° ; 

(b) \Jp^pDp{l) = A; 

(c) For all& y^ a (z P and all maps s : ct ^ [0,1], the intersection Si[cr, s) = Hpecr Dp[s[p)) is either 
empty or connected; 

(d) For all% ^ a <z P and all maps s : cr ^ [0, 1), the following implication holds: ®(cr,s) ^ and 
@{a,s') = % for alls' )- s implies that Qi{a,s)is reduced to a single point. 

Then, there exists a sequence of collapses which transforms Nrv{Dp{0) | p eP} intoNrv{Dp[l) \ p eP}. 

Proof. The nerve Kt = Nrv{Dp(f ) \ p ^P] changes while t increases continuously from to 1. Because 
of condition (a), the only changes that may occur are the disappearances of some simplices. Let Af 
be the set of simplices that disappear at time t. Suppose Af 7^ and let us prove that the deletion 
of simplices Af from Kt is a collapse for all t e [0, 1). Generically, we may assume that the set of 
simplices Af has a unique inclusion-minimal element CTmin. We will explain in the appendix how 
to get rid of this genericity assumption. Assuming we are in this generic situation, we do not need 
anymore conditions (c) and (d) but can replace them with the weaker conditions (c') and (d') obtained 
by considering constant maps for 5 and s'. Since f]p^^ Dp[t) 7^ and f]p^u . Dp[t + 17) = for all 
< J7 < 1 — f , condition (d') implies that f]p^u . Dp[t) = [a] for some a ^ A. It is easy to see that 
Af has a unique inclusion-maximal element auisx = {p ^ P \ a ^ Dp[t)}. Hence, Af consists of all 
cofaces of CTmin and these cofaces are faces of a max- To prove that removing Af from Kt is a collapse, 
it suffices to establish that amm 7^ CTmax- We proceed in two steps: 

Step 1: Let us prove that a lies on the boundary of Dp[t) for all p e CTmin- For this, define ctq = 
[p e CTuiin I a e dDp[t)}. We start by establishing that ctq 7^ 0. Indeed, suppose for a contradiction 
that a belongs to the interior of Dp{t) for all p e (J mm- This implies that for some 17 > 0, we have 
B[a,ri) c npeo-min ^P^^^ ~ ^^^' yielding a contradiction. Let us now prove that ctq = o-min- Suppose 
for a contradiction that ctq is a proper subset of CTmin- Consider ry > such that B{a,rj) c Dp{t) for all 
p e CTmin \ ctq. We have 

ae P Dp(f)n5(a,)7)c Q Dp{t) = {a\. 

and therefore f^p^^j Dp{t) is contained in {a} U {R'^ \ B{a, rj)) which consists of two connected com- 
ponents. The only possibility for (c') to hold is that f^^^ Dp{t) = {a}. Because of (a), for all t < t' <l, 
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we get that f] ^^^ Dp[t') c f] ^^ Dp[t)° = 0. It follows that ctq disappears at time t and the minimal- 
ity of amin implies that ctq = CTmin, yielding a contradiction. Thus, a lies on the boundary of Dp[t) for 

all p e CTmin. 

Step 2: Let us prove that CTmax 7^ CTmin- By condition (b), we have a E A = {Jp^pDpil) and therefore 
fl belongs to -0,^(1) for some q ^ P. Since f < 1, condition (a) implies that a E Dq{\) c Dq[t)° and 
therefore ^ e CTmax- On the other hand, a ^ dDq[t) and therefore (7 ^ CTmin- h follows that a^^^^ 7^ 

O-min- □ 

Remark. Somewhat surprisingly, condition (c) of Lemma 10 is weaker than the condition required 
by the Nerve Lemma for guaranteeing that the simplicial complex Kt = Nrv{Dp(r) | p e P} is homo- 
topy equivalent to A at some particular value of f E [0, 1]. In particular, if the Nerve Lemma applies 
at time f = 0, that is, if npeo- ^p ^^^ ^^ either empty or contractible for all 7^ cr c P and if furthermore 
the four conditions of Lemma 10 hold, then Kt wiU have the right homotopy type for aU f e [0, 1] . 

6.3 Final 

For proving Theorem 3, we build a 2-parameter family of compact sets ® = {Dp{t) \{p,t)^Px [0, 1]} 
in such a way that if we let K[t) = Nrv{Dp(0 | p e P}, then CechAiBa) = K{0) and /(Nrv'g') = K[l). 
We then prove that this family meets the hypotheses of Lemma 10, implying that Cechyi(i^ a) can be 
transformed into f[Nrv'ia) by a sequence of collapses obtained by increasing continuously t from 
to 1. To define the family ®, let us first associate to every point p e P the set 

fcienchers„(C,) iff-\p) = {v], 
^ \{p+,p-} iir\p) = &, 

where p+ and p~ are two points which are symmetric with respect to p and chosen such that B[p^, a)n 
B{p~,a) = 0. We then set 

Dp{t) = An P B{{l-t)p + tz,a). 
zeZip] 

Let us first check that CechA{P,a) = K{0) and /(Nrv-g") = K{1). We claim that Z{p) 7^ for all p e P. 
Let us consider two cases. First, if f'^ip) = 0, then by definition Z[p) = {p'*',p~} 7^ 0. Second, 
if f~^ip) = [v], then Z[p) contains at least p since Cp c B°[p,a) by hypothesis. Thus, Dp(0) = 
Anp^g2^pjB(p,a) = AnB{p,a) and ^(0) = Cech^(J!a). On the other hand, we have Dp{l) = An 
r\zez(p) ^^^' '^^ which we can rewrite as 

^ JAnHulUC) i{f-\p) = {v}, 

Thus, K[l) = /(NrvS^). By construction, each cell Dp[t) is the intersection of the shape with a collec- 
tion of balls whose center setZ(p) is either empty or compact and the collection of cells {Dp[t) | p e P} 
cover the shape. Thus, combining Lemma 7 and the Nerve Lemma, we get that NjrvK[t) ^ A for all 
t e [0, 1] . We are now ready to prove a stronger result, namely that as t increases continuously from 
to 1, the only changes that may occur in K[t) are collapses. 

Proof of Theorems. It suffices to establish that the family ® defined above satisfies conditions (a), 
(b), (c) and (d) of Lemma 10. 
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(a) Let us prove that for all < f < f ' < 1 and all peP,we have Dp{t') c Dp{tf. If /-i(p) = 0. this is 
easy to see. Suppose f~^{p) = {v}; see Figure 6, left. Since Cy c B°[p, a), we can find 5 > such that 
for every point x in B{p, 5), the inclusion Cy c B(x, a) holds. Equivalently, B[p, 5) c ClencherSa(Ci,). 
Consider the half-line Lu emanating from p in direction u. Because ClencherSa(Cy) is convex and 
compact, it intersects the half-line Lu in a closed line segment pzu with ||p — z^ || > 5. Since, for every 
point X on the closed line segment z-^uZu: we have B[Zu,a)n B[z-u,a) c B[x,a), it follows that one 
can rewrite Dp[t) as 

Dp{t) = An P Biil-t)p + tZu,a). 

Write p^ = il-t)p + tZu for short and observe that ||p^'-p^|| = (r'-r)||/?-z„|| >it'-t)5. Hence, as 
we move from z^u to Zu on the line segment z-f^z^, we meet the five points pi„, pL^, P, P^ and p^' 
in this order and the distance between the first two points and the last two points in the sequence is 
lower bounded by (t' — ?)5. We will now use the following observation. If we consider a point x on the 
line segmental suchthat ||a— x|| >e and ||fo — x|| >e, then B[a,a)nB[b,a)<z B{x, \J a^ — s^). Setting 
«o = V «^ -{f - 1)^5^, we thus get that B{p'_lu , a)nB{p'j^,a) c B[pLu, «o)n 5(p^ , ao) and since ao < a 

Dp{t') c ^n P B{{l-t)p + tZu,ao) c Dp{t)°. 



B°{p,a) 





Figure 6: Notation for the proof of Theorem 3. 

(b) Clearly Up6P-Dp(l) = A. 

(c) Given a czP and a map 5 : cr ^ [0, 1], we introduce the intersection of balls with radius a 

ncr,s]=f]Dpisip)) = Anf] f] Biil-sip))p + sip)z,a]. 

pea peazeZip) 

By Lemma 7, the intersection ®(cr, 5) is either empty or connected. 

(d) Consider a cz P and a map 5 : cr ^ [0, 1) such that ®(cr, 5) 7^ 0. Let us prove that if ®(cr, s') = for 
all maps s' :a ^> [0, 1] with 5' >- s, then ®(cr, s) is a singleton. Assume, by contradiction, that ®(cr, s) 
contains two points xi and X2 and let us prove that we can find a map s' : o" — > [0, 1] such that s' >- s 
and ®(cr,5') 7^0. Take a' such that a < a' < Reach (A). Since An Hulla'({xi,X2}) contains both xi and 
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X2, it is non-empty and therefore contractible by Lemma 7; see Figure 6, right. In particular, there is a 
path connecting the points xi andx2 in^nHuna'({xi,X2}). This path has to intersect the largest ball 
B contained in Hull„'({xi,X2}) and therefore A n B 7^ 0. For ^>0 sufficiently small we have 

AnB®'- c ^nHull„({xi,x2}) c ®(o-,5). 

By moving slightly the centers of the balls defining @(cr, 5), that is, by replacing the map s by a map 
s' >~ s such that s'[p) — s[p) is small enough for all p e cr, we get a new set &{(J, s') that still contains 
B. Since 7^ A n 5, we thus get that ®(cr, s') 7^ 0, reaching a contradiction. D 

7 Nicely triangulable spaces 

Given a space A and a sample P of A, we are seeking a sequence of collapses that transform the 
Cech complex of P with scale parameter a into a triangulation of A. We recall that a triangulation 
of A is a simplicial complex whose underlying space is homeomorphic to A. If A has a triangulation, 
then A is said to be triangulable. In particular, we know that smooth manifolds are triangulable [16] . 
Unfortunately, the proof involves barycentric subdivisions whose dual meshes are not likely to have 
convex cells and therefore have little chance of being a-robust coverings. And yet, we know that if 
a triangulation T of a space A is the nerve of an a-robust covering of A, then the previous section 
provides conditions under which Cech^(J)a) can be transformed into T by a sequence of collapses. 
This raises the question of whether, given a space A and a scale parameter a, it is possible to find a 
triangulation T of A which is the nerve of some a-robust covering of A. In this section, we focus on 
the question and exhibit spaces enjoying this property. 



^nHull„(C„) 




An[Conv(C„)]®''°'' 



HStTiv]) 



Figure 7: Left: the collection of disks and the collection of Voronoi regions both form an a-robust covering of the 
plane. Right: A triangulation is nice in our context when, among other things, it is the nerve of a collection of 
cells Ci, with size p such that An [Conv(C„)]®''°'' c ft(Str(i')) for some rjo > 0. This property will be preserved by C^' 
diffeomorphisms for p small enough. 

As a warm-up, we study the case A = R^; see Figure 7. Consider a Delaunay triangulation T of 
M^ with vertex set V and write C,, = {x e M^ | ||x — v\\ < \\x — u\\ for all u ^V] for the Voronoi cell of 
V ^V. Setting "iS = {Ci,\ v ^V} for the collection of Voronoi cells, we have that T = NrvS^. If all angles 
in T are acute, then the Voronoi cell Cy is contained in the star of v and so is Hulla(C[,) for a large 
enough. In particular, by choosing carefully V and a, we can ensure that "io is an a-robust covering of 
the plane. Writing Str(y) for the set of simplices of T incident to v and Conv{X) = Hull+c>o(^) for the 
convex hull of X, this motivates the following definition. 
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Definition 2 (nice triangulation). Let p and 5 be two positive real numbers. A triangulation T of 
AczR'^ is said to be [p , 5)-mce luith respect to [h,"^) ifh is a homeomorphism from \T\ toA, "^ = {Cv \ 
V &V} is a finite compact covering of A such that Nrv '^ = T and the following conditions hold: 

(i) h{a) c lJj;go. Ci, for all simplices cr e T; 

(ii) Cu<zB°ihiv),p)forallveV; 

(iii) Ar][Conv{Cv)f^ c hiStriv)) for all v ^ V . 

If ris a (p, 5)-nice triangulation oiA with respect to [h, "ia), then the next lemma gives conditions 
on a, p and 5 which guarantee that "ia is an a-robust covering of A. 

Lemma 11. Let A he a compact set ofW^ and suppose T is a{p ,5)-nice triangulation of A with respect 
to{h,'i§). Then'iS isana-robust covering of A whenever the following two conditions are fulfilled: (1) 
p <a and (2) a — y a^ — p^ <5. 

Proof Suppose '^ = {Cv\v^V} and let v^V. By Lemma 13, Hull„(C„) c [Conv(C„)]®^; see Figure 7, 
right. It follows that Cv c AnHulla(C(,) c h[StT{v)) from which we deduce the sequence of inclusions 

r = Nrv'^ c Nrv{AnHull„(Ci,) \veV](Z Nrv{h(Str(i^)) \veV} = T. 

The nerves on the left and on the right are equal, showing that Nrv'i^ = Nrv{AnHulla(Cy) \ v ^V}. D 

Observe that if T is a (p, riop)-nice triangulation oiA for some 170 > 0, then conditions (1) and (2) 
of Lemma 11 are satisfied for 5 = rjop as soon as p is small enough. Of course, the difficult question 
is whether such a triangulation T can always be found for arbitrarily small p . 

Definition 3 (nicely triangulable). We say that A c M^ is nicely triangulable if we can find Pq > and 
rio>0 such that for allO<p <po, there is a (p , rjop )-nice triangulation of A. 

Theorem 4. SupposeAczR'^ is nicely triangulable. ForeveryO<a<Reach[A), there exists sq > such 
that for all finite point set P czM.'^ and allO< e < sq satisfying A c P®'^, the complex CecliAiP, a) can be 
transformed into a triangulation of A by a sequence of collapses. 

Proof. By definition, we can find po> and ryo > such that for all < p < po, there is a (p, J70P)- 
nice triangulation T of A with respect to [h,'i§). Let us choose p small enough so that p < a and 
a— \] a? — p^ < rjop. Lemma 11 then implies that 'ig' is a compact a-robust covering of A. Set e{T,h) = 
|inf||/z(yi) — /i(i'2)|| where the infimum is over all pairs of vertices vi y^ V2 of T and let ^o be the 
minimum of e{T, h) and a — p. Consider a function / : Vert(r) — > P that maps each vertex i^ to a 
point ofP closest to h[v). Note that / is injective, \\h{v) - f[v)\\ <e and Cy cB°(/(i;),a) for all veV. 
Applying Theorem 3 yields the existence of a sequence of collapses from Cech^(i^a) to /(T). D 

We now establish that the property of being nicely triangulable is preserved by C^-^ diffeomor- 
phisms between manifolds. 

Theorem 5. LetM andM' be two compact C^-^ k-manifolds embedded respectively wR*^ and'R.'^' and 
^ : M ^> M' a C^-^ diffeomorphism. M is nicely triangulable if and only ifM' is nicely triangulable. 
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The proof is given in the Appendix. We conclude the paper with a theorem that provides a few 
templates of nicely triangulable manifolds. Thanks to Theorem 5 the property for a template to be 
nicely triangulable is inherited by all manifolds C^'^ diffeomorphic to it or equivalently by all mani- 
folds with a positive reach homeomorphic to it [14]. 

Theorem 6. The following embedded manifolds are nicely triangulable: 

1. The unit 2 -spheres'^ = {x = (xi,X2,X3)eIR^ | X!;=i-^f — ^l' 

2. The flat torus T^ = {x = (xi,X2,X3,X4)eIR4 | Xj +x| = 1 andx^ + x^ = 1}; 

3. The m -dimensional Euclidean space R™ , embedded in E*^ for some m <d. 

Proof. For A e {S^, T^, M*" }, we proceed as follows. We build a triangulation T parameterized by some 
integer n and consider a map h:\T\^>A. The integer n will control the size of elements in h[Ty. the 
larger n the smaller the image of simplices under h. We then consider the barycentric subdivision K 
of T and associate to each vertex y of T the cell Cy = IJo-gst (u) H^')- The collection of cells Cj, forms 
a covering "ia of A. In the three cases, it is not difficult to see that we can find J70 > such that T is 
(p, ?7op)-nice with respect to [h, "iS) for some p > 0. Furthermore, the value of p can be made as small 
as desired by increasing n. We thus conclude that A is nicely triangulable. Below, we just describe 
how T and h are chosen in each case. 

1. §2 is nicely triangulable. We start with an icosahedron centered at the origin and subdivide each 
triangular face into 4" equilateral triangles. The resulting triangulation T is then projected onto the 
sphere, using the projection map h : \T\ ^ E>^ defined by h[x) = j^tt. Let V be the vertex set of T. 
Notice that all vertices v e^V have degree 6 but the 12 vertices in the original icosahedron which have 
degree 5. As n increases, the star of each vertex v projects onto the sphere in a set that becomes 
arbitrarily close to a regular planar k-gon centered at v where k is the degree of v. 

2. T^ is nicely triangulable. The map // : M^ ^ fZ defined by //(s, t) = (cos 5, sins, cos t,sm t) is locally 
isometric and its restriction h : [0, 271)^ — > T^ is an homeomorphism. The idea is to build a periodic 
tiling of M^ made up of identical isosceles triangles as in Figure 8, left. Let a and b be the respective 
height and basis of the triangles. Consider two integers n and k such that na = kb = In. Taking 
fc = [^nj we get that the ratio ^ tends to ^ as n — > +00. Thus, the map H turns the periodic tiling 
of M? into a triangulation of T^ whose triangles become arbitrarily close to equilateral triangles with 
edge length b as n^> +00. 

3. M™ is nicely triangulable. We start with a cubical regular grid and define T as the barycentric 
subdivision of that grid; see Figure 8, right. Precisely, for each cell in the grid, we insert one vertex at 
its centroid. So each edge is subdivided into 2 edges sharing the inserted vertex. We then recursively 
subdivide the cells by ascending dimension. Each cubical k-cell has 2k cubical (fc — l)-cells on its 
boundary. We subdivide each fc-cell as a cone whose apex is the inserted vertex and whose basis is 
the subdivided boundary of that cell. We claim that all stars in T are convex. Indeed each vertex in 
T is the centroid of an initial cubical cell of dimension between and m. Consider the vertex v that 
was inserted at the center of the fc-dimensional cubical cell Dp and let us describe the set of vertices 
Vv in the link of v in T. The vertices of V^ can be partitioned in two subsets. The first subset contains 
vertices in the /c-flat that supports Du while the second subset contains vertices in the [d — fc)-fiat 
passing through v and orthogonal to the fc-flat supporting D^. The vertices in the first subset lie on 
the boundary of D^ and the vertices in the second subset lie on the boundary of a (m -k)-cube. Since 
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in both flats of respective dimension k and m — k the vertices in VJ, are in convex position, it results 
that vertices in Vy are in convex position in M™ . As a result, it can be proved (details are skipped) that 
the star of v is the convex hull of V^, . Finally, we let h be the identity map and n the inverse of the size 
of the grid. D 
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Figure 8: Triangulating T^ (left) and R^ (right). 



8 Conclusion 

In this paper, we introduce the notion of nicely triangulable shapes. We prove that the Cech complex 
of a point set at scale a can be simplified by collapses into a triangulation of a nicely triangulable 
shape, assuming the point set samples well enough the shape and the scale parameter a has been 
well chosen. Precisely, the sampling conditions and values of a are exactly the same as the ones 
required by the reconstruction theorem of Niyogi, Smale and Weinberger [15] for guaranteeing that 
the a-offset of a point set deformation retracts to the shape that the points sample. We conclude with 
a few open questions: 

(1) Our result assumes the shape to be nicely triangulated. In this paper, we list a few simple spaces 
which enjoy this property. The list includes all surfaces C^-^ diffeomorphic to a sphere. Is it pos- 
sible to extend this list to a larger class of manifolds? In particular, we conjecture that surfaces 
embedded in R^ with a positive reach are nicely triangulable. 

(2) Unfortunately, our proof is not constructive. Indeed, the order in which to collapse faces is de- 
termined by sweeping space with a f -offset of the shape for decreasing values of t. Since the 
common setting consists in describing the shape through a finite sample, the knowledge of the 
f -offsets of the shape is lost. Nonetheless, is it possible to turn our proof into an algorithm? 
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A Two geometric lemmas 

The following lemma is useful for establishing Lemma 3. 

Lemma 12. Let X c M^ be a non-empty compact set and B[c,p) its smallest enclosing ball. For all 
points z ^R'^ and all real numbers a>p, the following implications hold 



(i) X c B(z, a) =^ B(c, a - y/ a^ - p'^) c B{z,a); 



(ii) B(c,a-^a2-p2)cB°(z,a) =^ XnB°(z,a)7^0. 

Proof. To establish (i), assume for a contradiction that B(z, a) does not contain B(c, a — y a^ — p^) or 

equivalentlythat ||c — z|| > \] o? —p'^. This implies that the smallest ball enclosing B(c,p)nB(z, a) has 
radius less than p. Since this intersection contains X, this would contradict the definition of p as the 
radius of the smallest ball enclosing X. 

It is not hard to check (by contradiction) that the center of the smallest ball enclosing X lies on 
the convex hull of X. It follows that for any half-space H whose boundary passes through c, the 
intersection Xn i?n B[c,p) is non-empty. To establish (ii), we may assume that c 7^ z for otherwise 
the result is clear. Let H be the half-space containing z whose boundary passes through c and is 
orthogonal to the segment cz; see Figure 9, left. If B(c,a— \/a^ — p^) c B°(z,a), then Hn B{c,p) c 
B°(z, a). Since Xr\Hr\ B[c,p)j^ 0, it follows that Xn B°(z, a) 5^ 0. D 

The following lemma is useful for establishing Lemma 11. 

Lemma 13. LetXd M^ be a non-empty compact set and B[c,p) its smallest enclosing ball. For all a 
and 5 such thata> p and a — y a^ — p^ <5, the following inclusion holds: Hulla{X) c [Conv(X)]® . 

Proof. See Figure 9, right. Consider a unit vector u e S'^~^ and let Lu be the half-line emanating from 
c in direction u . Let B^ denote the ball with radius a centered on L„ containing X and whose center 
is furthest away from c. By construction, X is contained in the intersection of the two balls B^ and 
B{c,p). The boundary of B^ n B{c,p) consists of two spherical caps and we let C^ be the one lying 
of the sphere bounding B^ . Observe that X has a non-empty intersection with C^ and for all p >a, 
the ball B^ intersects C^. The largest distance between a point of C'^ and B^ is upper bounded by 

the height of C" which is less than of equal to a - \/a^-p^ < 5. We thus get that B^ c [B^]®^. 
Considering this inclusion over all directions u for p = -l-oo yields the result. D 



B C ^ ' ^ diff eomorphism preserves nicely triangulable manifolds 

The goal of this section is to prove Theorem 5. Let us provide some additional definitions. A C^-^ 
function is a differentiable function with a Lipschitz derivative. A C^'^ structure on a manifold is 
an equivalence class of atlases whose transition functions are C^-^ Finally, C^'^ diffeormorphisms 
betwen C^'^ manifolds are defined accordingly. In Theorem 5, we restrict our attention to shapes 
which are C^-^ compact manifolds without boundary embedded in M^ and whose embedding are 
themselves C^-^ (for the natural smooth structure of M''). We will say that such shapes are compact 
C^'^ manifolds embedded in R'' for short. A compact manifold embedded in R*^ is C^'^ if and only if 
it has a positive reach [14] . 
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Figure 9: Notation for the proof of Lemma 12 (left) when ||c-z|| = \l a^ - p'^ and for the proof of Lemma 13 (right). 

Proof of Theorem 5. Let x E M. Since M is a compact C^-^ fc-manifold embedded in W^, there exists 
a A; -dimensional affine space Tm{x) c M'^ tangent to M at x. Let tt^ : M — > TmIx) be tiie orthogonal 
projection onto the tangent space Tm{x) and let 7r^,_^, : M' — > TM'i^ix)) the orthogonal projection onto 
Tm'{^{x)). Since M and M' are compact, we can find two constants K and K' independent of x such 
that: 

VyeM, ||y-7r,(y)||<^||y-x||2 (3) 

Vj/eM', Il3/-<(,)(>')ll<^'lly-$(x)||2 (4) 

Given to > 0, we consider the open set Ux = Mn B°(x, to) and adjust to in such a way that 

1. The restriction nx:Ux^>7ixiUx) is an homeomorphism for all x e M; 

2. The restriction tt^,^, : $([/x) — > 7r^,^J$([/x)) is also an homeomorphism for all x e M. 

The collection of pairs {{Ux,7ix)}xeM forms an atlas in the C^'^ structure of M. Similarly, the collection 
of pairs {($(C/x),7iL^0}xeM forms an atlas in the C^'^ structure of M'. Let 5^m(x) be the linear space 
associated to TmIx) and denote by D$x the derivative of $ at x, seen as a linear map between ■9'm{.x) 
and 5^M'(^(^))- Since M is compact, there is a constant Kq, independent of x such that: 



Vy e M, my) - $(x) - D<^x[7ix[y) - x)\ < ^$||y - x||2. 



(5) 



and two positive numbers kz > ki > 0, again independent of x by compactness of M, such that 

Vwe57k(x), K1||U||<||D$:,(M)||<R-2||U||. (6) 

Consider the affine function $x : Tm(-^)^ rM'($(x)) defined by $x(y) = $(-^) + £^$x(y—^)- Combining 
Equations (3) (4) and (5), we can find a constant L$ independent of x such that for all t < to and all 
compact sets A c Mn B(x, t) : 



dni^x ° 7lx{A),7lL.o^{A)) < Lq,t 



(7) 



Now, assume that M is nicely triangulable and let us prove that M' is also nicely triangulable. 
By definition, we can find po > and rjo> Q such that, for all < p < po, there is a (p, }7op)-nice 
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triangulation T of M with respect to some [h, '^). Suppose "iS = {Ci,\ v ^V] and consider tiie covering 
S^' = {$(C[;) I V e V}, tiie iiomeomorpiiism h' = ^ o h : \T\ ^> M' , the real numbers p' = 2k2P and 
^70 ~ '^"^°2^. — ^- Let us prove that by choosing p small enough, T is a [p', rJQp'ymce triangulation of 
M' with respect io{'€' ,h').ln other words, we need to check that conditions (ii) and (iii) of Definition 
2 are satisfied for "^ = '^', h = h' , p = p' and 5 = rj^p'. Take v ^ V and set x = h[u), C = Cy, 
S = h{StT{v)). 

(ii) By definition of T, we have x e C c B°(x,p). Taking the image of this relation under $ and 
choosing p > small enough, we get that $(x) e $(C) c $(B°(x,p)) c B°($(x),p'). The last inclusion 
is obtained by combining Equations (3), (5) and (6). 

(iii) Let us choose a positive real number p < min{po, y} small enough to ensure that r]Q>0 and 
let us prove that M' n [Conv(C)]®''°'^ c S. By choice of T as a [p,riop)-mce triangulation of M with 
respect to (h,'^), we have that M n Conv(C)®'"'/' c S. Furthermore, C c B{x,p) and S c B[x,2p). 
Thus, by choosing p < j, we have SczUx and 

[/x n Conv(C)®'"'^ c S. 

Taking the image by the homeomorphism tIx : Ux ^> '^xiUx) on both sides and using UxiAn B) = 
7ix[A)nnx[B)weget 

7lx{Corw{Cf'^oP)^7Zx{S). 

Let 5fc(0. ?") denote the fc-dimensional ball of ■:^mM centered at the origin with radius r. Writing 
A® B = {a + b \ a ^A,b ^ B} ior the Minkowski sum of yl and B, it is not too difficult to prove that 

TixiA®^) = TixiA) ® 5fc(0, 5). It follows that 

TixiConviC)) ® BkiO, r]op) c nx{S). 
Taking the image under ^x on both sides we get 

$x ° 7rx(Conv(C)) ® D^xBidO, rjop) c $x ° nx[S)). 

Let 5^(0, r) denote the fc-dimensional ball of .5^m'(^(^)) centered at the origin with radius r. Using 
Equation (6) we get that B'^{Q,KiirjQp) c D^xBkiO.rjop). Since $x and tIx are both affine, so is the 
composition and therefore ^x ° 7rx(Conv(C)) = Conv($x ° TixiC)). It follows that 

Conv($x ° 7ix{C)) ® B[{0, Kir]op) c $x ° ^x{S). 

Recalling that C c B[x,p) and S c B(x,2p) and combining the above inclusion with Equation (7) we 
obtain 

Conv(7r; o $(C)) ® Bfc(0, KH70P - 5Up^) czn'^o $(S). 

Interchanging Conv and n'^, noting that rj^p' = Kirjop — 5L$p2 and using 7i'^[A®^) = 7r^(A) ® Bjt(0, 5) 
we get 

7r;(Conv($(C))®''>') c tt^ o $(S). 

Since tt'^ : $(C/x) — » ^$(x)(^^(^^-*:)) i^ homeomorphic, we thus obtain M' n Conv(${C))®''(i^ c $[S) as 
desired. D 
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C Getting rid of genericity assumptions 

In the proof of Lemma 10, we consider a 1-parameter family of nerves, Kt = Nrv{Dp[t) | p e P| for 
t e [0, 1]. Letting Af be the set of simplices that disappear at time t, we assume that the following 
generic condition is satisfied: 

(•) If Af 7^ for some t E [0, 1), then Af has a unique inclusion-minimal element CTmin- 

If not the case, the idea is to apply a small perturbation to the family ® = {Dp[t) | p e P| which 
will leave unchanged Ko and Ki and such that after perturbation (1) ® will still satisfy the hypotheses 
of the Lemma 10; (2) the generic condition {*) will hold. 

Definition 4. Two simplices a i andcTz are said to be in conjunction at time t if they are both inclusion- 
minimal elements of At for some t e [a,b). 

Consider two simplices that are in conjunction at time t, say cri and ctz- Suppose p ^ a\ and 
p ^ (72. Consider an increasing bijection ij) : [0, 1] -^ [0, 1] that differs from identity only in a small 
neighborhood of t that does not include any other event times. By replacing Dpby DpOip, we change 
the time at which cri disappears while keeping unchanged the time at which ctz disappears. After this 
operation, Ui and az are not in conjunction anymore. Furthermore, the operation does not create 
any new pair of simplices in conjunction. By repeating this operation a finite number of times, we 
thus get a new collection as required. 

The generic condition in the proof of Theorem 2 is removed in a similar way. 
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